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Dark matter (DM) particle annihilation or decay can produce monochromatic 7-rays readily 
distinguishable from astrophysical sources. 7-ray line limits from 30 GeV to 200 GeV obtained 
from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a 
selection based on requirements for a 7-ray line analysis, and integrated over most of the sky. We 
obtain 7-ray line flux upper limits in the range 0.6 — 4.5 x 10 _9 cm _2 s _1 , and give corresponding DM 
annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed. 

PACS numbers: 95.35.+d, 95.85.Pw, 98.70.Vc 



Introduction: Annihilation or decay of dark matter 
(DM) particles can give rise to Standard Model parti- 
cles, including 7-rays. The predicted 7-ray flux is usu- 
ally weak compared to astrophysical sources, and strik- 
ing experimental signatures are valuable in recognizing 
a signal. A distinctive signal would be the detection of 
a narrow 7-ray line [Il-Q originating from dark matter 
particles annihilating or decaying into 7X, where X can 
be another photon, a Z-boson, a Higgs-boson, or a non- 
standard Model particle. A 7-ray line search using data 
from EGRET in the energy range 1 — 10 GeV has been 
presented in [8] . We present here for the first time results 
from 11 months of Fermi Large Area Telescope (LAT) 
data in the energy range 20-300 GeV and obtain limits 
for 30-200 GeV. 



DM particles of mass m x annihilating into produce 

monochromatic 7-rays of energy E 1 = m x (l — j^^J, 

where mx is the mass of X. For decays x ~ ► lX, the 
equation gives the photon energy after the substitution 
m x — > m x /2. Since dark matter is strongly constrained 
to be electrically neutral, it has no direct coupling to pho- 
tons. The process x{x) ~~ ^ jX thus occurs only at higher 
orders, and with a branching fraction that is typically 
only 10~ 4 — 10 _1 compared to the total annihilation or 
decay rate. However, a strong line feature, or a strong 
line-like feature such as a sharp cut-off at a particular 
energy, can be expected in models with supersymmetry 
0-13 1, an extended Higgs sector 0], from final state ra- 
diation IB" 12 1 , from gravitino decay 18], or in scenarios 
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with non-thermal WIMP production 

LAT Data Selection and Analysis: The LAT is a pair- 
conversion telescope that combines silicon-strip/tungsten 
trackers and hodoscopic CsI(Tl) calorimeters into a 4x4 
array of 16 identical modules. The tracker is covered by 
a segmented anti-coincidence detector (ACD). Including 
the tracker, the LAT presents 10 radiation lengths for 
normal incidence. The depth, segmentation, and wide 
field of view of the LAT enable its high-energy reach. 
Separation of the very large charged cosmic-ray back- 
ground from 7-rays is achieved via the combination of 
the data acquisition trigger, on-board event software fil- 
ter system, and extensive ground processing of the data. 
Details of the LAT, and the data analysis are given in 
1201 . An account of how the LAT is calibrated on orbit 



is presented in [21] . The LAT nominally operates in a 
scanning mode that covers the whole sky every two or- 
bits (~ 3 hours). The analysis described here uses data 
taken in this scanning mode from Aug 7, 2008, to July 21, 
2009, corresponding to an average exposure of 3.3 x 10 10 
cm 2 s. 

Events are selected for this analysis only if they pass 
additional cuts to the Pass6V3 diffuse class cuts, i.e., the 
cleanest photon sample in the currently public Fermi data 
release 2(J 22 j. These additional cuts are: a) a small av- 
erage charge deposited in the tracker planes (veto against 
heavy ions); b) a transverse shower size in the calorimeter 
within a size range expected for electromagnetic show- 
ers (veto against hadronic showers and minimum ioniz- 
ing particles). Cuts a) and b) dramatically reduce the 
charged particle background at the loss of some effective 
area, yielding a 7-ray efficiency > 90% relative to that 
of Pass6V3 diffuse class. These cuts remove charged par- 
ticle backgrounds that would worsen our upper limits, 
and that can even yield structures that might be inter- 
preted as 7-ray lines. These cuts are a subset of those 
used in the LAT Measurements of the Isotropic Diffuse 
Gamma- Ray Emission [23] . 

In addition, we use only one of the three energy mea- 
surement methods from the LAT standard analysis. The 
method used is the LAT profile method [20| . in order 
to not introduce energy dependent artifacts that arise 
from switching between methods over the energy range 
discussed here. In the profile method, the energy is ob- 
tained from a fit to the longitudinal shower profile while 
considering the transverse profile. The exclusive use of 
the profile energy method led to no additional reduction 
in efficiency. The instrument response functions (IRFs) 
are recalculated for this data selection and for the use of 
the profile method. 

The resulting energy resolution averaged over the LAT 
acceptance is 11% FWHM for 20-100 GeV, increasing to 
13% FWHM for 150-200 GeV. The photon angular reso- 
lution is less than 0.1° over the energy range of interest 
(68% containment). The absolute calibration and energy 
resolution of the LAT was determined by comparing with 



e~ beam test data, taken at CERN in a secondary e~ 
beam with energies up to 300 GeV using a special cali- 
bration unit made u p of flight spare towers and ACD tiles 
(not the LAT itself ) [24|. The energy resolution measured 
in the beam tests is in agreement with the predictions 
from the Monte Carlo simulator based on GEANT4 that 
was used to define the IRFs (GLEAM [H). Also, the 
systematic error on the absolute energy of the LAT was 
determined to be -10 + 5% for 20-300 GeV. 

The systematic uncertainties for the exposure over this 
energy range are ±20% based on the extrapolation of 
studies comparing the efficiencies of analysis cuts for data 
and simulation of observations of Vela [25|. The Vela 
studies cut off at 10 GeV, and measurements of the sys- 
tematic errors above 10 GeV will be made when sufficient 
statistics are available from high-energy sources. Thus, 
the exposure systematic errors above 10 GeV that we 
quote for this study have not been fully validated. We 
believe that any reasonably projected uncertainty would 
not have a significant impact on the interpretation of the 
limits presented here. 

In our search for lines we use a region of interest (ROI) 
that covers most of the sky: an all-sky ROI with the 
Galactic plane removed (i.e., \b\ > 10° as the Galactic 
plane is very bright in photons from gas interactions) 
plus a 20° x 20° square centered on the Galactic center 
(GC) and aligned on the (£, b) grid of the Galactic coor- 
dinate system (Galactic coordinates in degrees are used 
in this paper). Though increasing the photon flux lim- 
its averaged over the reported energy range by less than 
10%, including the GC gives significantly better theoret- 
ical line limits. For the highly point source rich region 
within 1° of the GC, no point source removal was done 
as this would have removed the GC. For the remaining 
part of the ROI, point sources obtained from the year-1 
catalog under development by the LAT team [26| were 
masked from the analysis using a circle of radius 0.2° 
centered on the measured point source position (conser- 
vative, considering the angular resolution above 20 GeV). 
This last cut removes about 0.4% of the solid angle and 
about 5% of the total photons. 

In searching for deviations from a locally-determined 
background, we use a sliding energy window with the 
window size adjusted to reflect the energy resolution. 
This results in limits for 30 < E 1 < 200 GeV us- 
ing data in the range 20 < E 1 < 300 GeV. The re- 
sponse of the LAT to a line feature in energy were de- 
termined from full detector simulations (GLEAM [20j|). 
Fig. [T] shows a binned representation of the fit and a 
close-up of the LAT line response function at 40 GeV. 
This fit also shows the largest line "signal" that was ob- 
tained in the reported energy region. We use an un- 
binned likelihood method, with the likelihood function 

ntot 

L{E\f,T) = n / ' S(Ei) + (1- /) • B(E i: T). Here, 
Ei denotes the measured energy of the i th photon. The 
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FIG. 1: A binned representation of a typical fit (done un- 
binned), here centered at 40 GeV, used to extract the flux 
upper limits presented in the tables. The fitting process is 
described in the text. In the main part of the figure, the 
lower (upper) dotted line is the signal (background) from the 
fit and the red(or black) line is the total fit. This fit also shows 
the largest line "signal" that was obtained in the reported en- 
ergy region. The inset shows a blow-up of the signal, which is 
the line energy response function, S(E), used in this fit, and 
is typical of line shapes for 20-300 GeV. 



parameters / and V are free and represent the signal frac- 
tion and the index of the power-law function, B (Ei, T), 
used to model the background. We require / > in the 
fit. The function S (Ei) models the signal shape, i.e., the 
LAT response for a line feature in energy averaged over 
the acceptance of the LAT. The confidence intervals are 
determined using the profile likelihood method (MINOS 
within MINUIT)[27], which provides two sided confi- 
dence intervals. The properties of this statistical method 
(coverage and power) have been thoroughly tested. At 
100 GeV, for example, the coverage is close to nominal 
for a range of true signal fractions from to 50%, and 
the power reaches 100% for signal fractions of about 10%. 
The method overcovers slightly due to the physical con- 
straint on the signal fraction, /. 

Results & Discussion: Table U shows flux limits as a 
function of photon energy that can be translated into 
bounds on the annihilation cross-section or decay life- 
time assuming a specific halo dark matter density pro- 
file, p(r). The monochromatic gamma-ray flux from dark 
matter annihilating into -fX with a cross-section (eru) is 



8-7T TT& 

otherwise. Here, 



£, where _/V 7 = 2 for X = 7 and N 7 = 1 



C 



db d£ ds cosbp 2 (r), 



(1) 



R% 



where the integral is over the ROI, r — (s 2 
2sRq cos £ cos 6) 1 / 2 , and Rp) — 8.5 kpc is the distance 
from the sun to the GC [28j. For decays, the flux is 
given by substituting in the equation for (av)/2m 2 , — > 



l/rm x , where r is the DM lifetime, and p 2 — > p in Eq. ([T]) 
for C. 

We consider three theoretically-motivated halo pro- 
files: the NFW profile, p NFW (r) = p s /[(r/r s )(l + r/r s ) 2 } 
with r s = 20 kpc [2{|, the Einasto profile, PEinastoM = 
p s expj — (2/a)\(r/r s ) a — 1]} with r s = 20 kpc and a = 
0.17 [30j, |31|, and the very shallow isothermal profile 
Pisothcrmai(» = p s /(l + (r / V s ) 2 ) with r s = 5 kpc [32|. We 
determine p s using p(Rq) — 0.4 GeV cm -3 [33j. Taking 
the mass of the Milky- Way halo to be - 1.2 x 10 12 M Q 
(see e.g. 34, l35j|). we determine maximum values for r 
of ~ 150 kpc for the Einasto and NFW profiles, and 
~ 100 kpc for the isothermal profile. 

Table |T] shows the cross-section and lifetime limits for 
the above profiles. We verified that there is only a minor 
dependence of the flux upper limits when changing the 
lower bound of \b\ in the range 8° < |6| < 15°, leaving the 
GC region fixed. The cross-section limits are sensitive to 
the halo profile. For the ROI used here, DM annihilation 
has a value of £ (Eq. J]}) for the Einasto profile that is 
40% larger than for the NFW profile, while for DM decays 
the C values are almost the same. This sensitivity is 
greater for cuspier profiles than those discussed here. For 
example the Moore density profile (3(| gives a factor of 
~ 6 stronger limits than the Einasto profile, using lower 
bounds on i and b for the Moore profile integration that 
correspond to a distance of 10~ 3 pc from the Galactic 
center. 

The limits on (crw) 77 ((erv)~ f z) shown in Table U are 
about one or more orders of magnitude weaker than 
the cross-sections expected for a typical thermal WIMP. 
However, there are several models in the literature that 
predict larger cross-sections and are constrained by these 
results. A WIMP produced non-thermally may have a 
much larger annihilation cross-section than a thermally 
produced WIMP and still produce the required DM relic 
density. An example is the "Wino LSP" model [19( that 
explains the recent positron measurement by PAMELA 



3Jj, and predicts (av)~ lZ — 1.4 x 10 ~ 6 cm 3 s 



at 



E. 



j ~ 170 GeV. Our results disfavor this model by about 
a factor of ~ 2 — 5, depending on the dark matter halo 
profile (see Table H]). Other models that are partially 
constrained include [38[ , while models that are only con- 
strained assuming a much cuspier profile include [39|, H(| ■ 
Dark matter decays into jX are ruled out for lifetimes 
below ~ 10 29 s, a limit that is largely independent of the 
dark matter halo profile. This constrains, for example, 
a subset of the lifetime range of interest for gravitinos 
decaying into mono-energetic photons [l8j |. 
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TABLE I: Flux, annihilation cross-section upper limits, and decay lifetime lower limits: 7-ray energies measured and cor- 
responding 95% c.l. upper limits (CLUL) on fluxes, for |6| > 10° plus a 20° x 20° square around the Galactic center. For 
each energy and flux limit, (cjv) 11 and (av) 1 z upper limits, and r 77 and r 7 z lower limits are given for three Galactic dark 
matter distributions (see text). The systematic error in the absolute energy of the LAT discussed in the text propagates to a 
— 20% + 10% systematic error on (eru) 77 , while for the decay lower limits the systematic error in the absolute energy of the 
LAT discussed in the text propagates to a 4-10% — 5% systematic error on r 77 . 
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